Preprint typeset in JHEP style - HYPER VERSION 



Inflationary Potential Reconstruction for a 
WMAP Running Power Spectrum 



James M. Cline, Loison Hoi 

Physics Department, McGill University 
3600 University Street, Montreal 
Quebec, Canada, H3A 2T8 

Abstract: The first year WMAP measurement of the CMB temperature anisotropy 
is intriguingly consistent with a larger running of the inflationary scalar spectral 
index than would be expected for single-field inflation. We revisit the issue of a large 
running spectral index, first by reexamining the evidence from the data, and then 
by reconstructing the inflationary potential, using an improved method based upon 
the Hamilton- Jacobi formulation. We note that a spectrum which runs only over 1.5 
decades of k space provides as good a fit to the CMB data as one which runs at all 
k, that significant evidence for running comes from multipoles / near 40, and that 
large running gives a better fit than a flat spectrum primarily if the tensor-to-scalar 
ratio r is large, r ~ 0.5, and the field values are at the Planck scale. This allows 
one to break the large degeneracy of potentials which would be consistent with the 
scalar power alone. Large running, should it be confirmed, is thus linked to a high 
scale of inflation and the possibility of seeing effects of tensor modes in the CMB and 
Planck-scale physics. Nevertheless, we show that the reconstructed inflaton potential 
is well-described by a renormalizable potential whose quantum corrections are under 
control despite the large field values. 
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1. Introduction 

In 2003 the Wilkinson Microwave Anisotropy Probe (WMAP) released the results 
of its first year of high-precision measurements of the temperature anistropy of the 
Cosmic Microwave Background (CMB) [p]]. While several anomalies of marginal 
statistical significance have been noted, one of the most visually striking is the trend 
of the spectral index of the scalar power n s , to run toward smaller values , shown in 
figure [I]. There it was determined that dn s /d In k ~ —0.1 albeit with large error bars. 
Like the other anomalies, this one also remains to be proven statistically significant, 
since a flat spectrum is also consistent with the data. Nevertheless, since the present 
data have opened the door to the possibility of large running, it is interesting to 
explore the theoretical consequences while we continue to wait for the next WMAP 
data release. 

It is somewhat surprising that the issue of running has not attracted as much 
attention as the low quadrupole since it is arguably as significant. Nevertheless 
there have been a number of papers which have treated running, starting even before 
WMAP @ . Refs. [|], H were the first to examine the likelihood of a running spectrum 
using modern CMB and Lyman a data. Subsequently to WMAP, ref. || found 
that the marginal evidence for running disappeared if one omitted the first three 
multipoles, a conclusion which we shall challenge below. Including data from the 
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Figure 1: Experimental limit on the 
scale-dependent spectral index, repro- 
duced from figure 2 of Q. Shaded areas 
and dashed lines indicate 68% and 95% 
confidence regions, respectively. 
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Figure 2: Probability distribution func- 
tion for running from combined WMAP 
and SDSS data, from Q. 



Sloan Digital Sky Survey galaxy survey slightly narrows the errors on dn s /d\nk, 
but not the central value, as shown in figures |1]-|2|. 

From the theoretical perspective, the —0.1 level of running is uncomfortably 
large, since dn s /d\nk is a second order effect in the slow-roll expansion; generically 
one expects dn s / dink < 1/ iV| with iV e ~ 60 e-foldings of inflation since horizon cross- 
ing. Ref. H considered which general classes of models could give large running, and 



methods of constructing potentials with large running. Refs. Q-JTlJ showed that 
there is a new consistency condition between the scalar and tensor power spectra in 
the presence of running. Ref. |T2j analyzed generically large running in conjunction 
with multifield inflation, and sudden changes of direction in field space. Supergrav- 
ity models with two-stages of inflation have been shown to provide large running 
. Classes of well-motivated particle physics models consistent with large running 



were presented in [fUjl ; inflation in noncommutative spacetime has been found to 



give significant running [|15| , |16 |, as well as inflation with a cutoff on transPlanckian 
fluctuations [ I7| . In |18| it was noted that although brane-antibrane inflation models 
predict smaller levels of running, dn s /dlnk ~ —0.01, nevertheless the correlation of 
dn s /dlnk with n s can provide an important test of the model if future experiments 
are able to reliably measure running at this low level. 

One possible reason for a decline of interest in running of n s is that recent Lyman- 
a forest data have moved the central value much closer to zero, and narrowed the 
confidence intervals. Ref. |1| finds that dn s /dhik = —0.003 ± 0.010 (see also p0| ). 
This determination would seem to put the question of large running to rest. How- 
ever, given the reasonable question as to whether systematic errors are completely 
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under control in the Lyman-a studies, we suspend judgment until there is further 
confirmation from WMAP or other CMB experiments. Moreover, Lyman-a data are 
sensitive to the smallest scales. It is perfectly possible that the power spectrum runs 
significantly at larger scales, but flattens out at smaller scales. Indeed, this is almost 
a necessity since it is impossible to sustain inflation over 60 e-foldings if the running 
remains large throughout, as we will discuss. 

There is a broader issue connected with increasingly precise experimental de- 
terminations of the inflationary power spectrum, namely, what is the shape of the 



inflaton potential which can generate a given spectrum |21], [22j? One of the points 
of this paper is that an expansion of the potential around a point using slow-roll 
parameters, and integrating the slow-roll approximation of the inflaton equations 
of motion, are inadequate for accurately reproducing the desired spectral shape, if 
it runs at the level which we are considering, or otherwise departs from the slow- 
roll approximation. In this work we present a simpler technique for reconstructing 
the potential, based on the Hamilton- Jacobi formulation of the inflaton equation of 
motion f2|, ||§. 

We begin in section |2| with a description of the running spectrum and its like- 
lihood, with our own analysis of CMB and large scale structure (LSS) data sets. 
We argue that the possibility of large running, though unproven, is still plausible, if 
applied to a restricted range of the power spectrum. In section |] we describe our new 
method of reconstructing the potential, based on the Hamilton- Jacobi method. A 
one-parameter family of potentials corresponding to the running scalar power spec- 
trum is constructed. In section [T3] we show how higher order corrections in the 
slow-roll expansion for the power spectrum can be efficiently incorporated. In sec- 
tion |3]4] we note that a given potential yields the desired power spectrum only for a 
specific initial velocity for the inflaton, which does not generally correspond to the 
slow-roll attractor value. In section [| we break the degeneracy between potentials 
found in the previous section by assuming a large tensor-to-scalar ratio, as is favored 
by the best fit parameters found in section ||. We discuss several aspects of this 
potential, including the duration of inflation which can be supported during the pe- 
riod of large running, and the extent to which the potential can be approximated by 
renormalizable terms in a quantum field theory. We give conclusions in section |5|. 



2. Running Power Spectrum 

The power spectrum with running spectral index is parametrized as 
P{k) = A exp 

The running spectrum which best fits the first year WMAP data was given in 
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at a pivot point of ko = 0.002 Mpc -1 , with a likelihood given by —2 InC = x 2 = 1426 
for 1340 degrees of freedom (number of data points minus number of parameters). 
The other cosmological parameters are listed in table |l|, where we show the best-fit 
running and nonrunning models which we obtain using the CosmoMC Markov Chain 
Monte Carlo program [5S| (see http://cosmologist.info/cosmomc). The \ 2 = 1426 
value for the running spectrum includes a penalty of A% 2 = 1.4 from the high value 
of the Hubble parameter, h = 0.814, relative to the Hubble Space Telescope value 
h = 0.72 [EG]. 1 All models shown are flat: = 1 — Q m . We note that a sizeable 
tensor-to-scalar ratio r (as well as a large optical depth r) is required to get the best 
fit with large running. If we set r = while keeping n s and dn s /d\nk fixed (but 
allowing other cosmological parameters to vary) the likelihood decreases, x 2 = 1450. 
There is a degeneracy between running and tensors; if we set r = and allow all 
other parameters to vary, including n s and dn s /dlnk, the best fit has x 2 — 1429, but 
smaller running, dn s /d\nk ~ —0.009. Although the statistics are not yet significant, 
there is the hint that large running and a large tensor ratio fit the data better. 
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Table 1: Our best fits for the running and nonrunning models using CosmoMC and the 
first year WMAP data. 



We have repeated the comparison between running and nonrunning models using 
data from other CMB experiments (ACBAR |27|, VSA @, CBI @), and galaxy 
surveys (2dF f30fl , SDSS @). (In all the data sets we also include the Hubble Space 
Telescope determination of the Hubble parameter, which contributes 156.25 x {h — 
0.72) 2 to the x 2 i leading to A% 2 ~ 1.4 between the running and nonrunning models 
using the WMAP data.) The best fit parameters for different combinations of these 
data are shown in table |^. Inclusion of other CMB data strengthen the evidence 
for running, while including large scale structure (LSS) data weaken the evidence. 
This is more evident with the inclusion of Sloan Digital Sky Survey (SDSS) data 
than with the 2dF data. 2 This is consistent with the evidence for running coming 
principally from the low-k part of the power spectrum, to which the CMB is more 



1 Nevertheless our value is somewhat smaller than the value obtained by WMAP for the running 
model §, x 2 = 1428. 

2 We have observed that there is a tendency for CosmoMC to get stuck in local minima of x 2 
when searching for the best fit parameters, especially when SDSS is included, so it is possible that 
we have not found the global minimum x 2 ■ 



- 4 - 



sensitive. If the effect turns out to be real, then the best fit to the data including 
large scale structure will require a spectrum which has some different behavior from 
eq. (|2.1j) at large k, as we will discuss below. 
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Table 2: Best fit parameters for both running and nonrunning spectral index models from 
CosmoMC. The first row is the result from the first year WMAP data alone; others are 
from the combinations of the WMAP data and the data sets specified at the first column, 
where CMB=ACBAR+CBI+VSA, mpk=2dF+SDSS. We used the July 2005 version of 
CosmoMC, and its default data sets. They are: the first year WMAP data, 2 < /tt < 900, 



2 < Zte < 450 (1348 data points); July 2002 ACBAR data [27|, 900 < I < 1950 (7 data 



points); 2000 and 2001 CBI data @, 700 < I < 1760 (8 data points); Feb 2004 VSA data 
|2§, 640 < I < 1700 (8 data points); June 2002 2dF data @, 32 data points; Oct 2003 
SDSS data 0, 0.016 < fc(h/Mpc) < 0.205 (19 data points). 



The relevant part of the spectrum is shown in figure [| the range of k values 
here is that which is actually called by the CAMB program (see http://camb.info/) 
in computing the multipole moments of the temperature anisotropy and polarization 
of the WMAP data. The scales k = 10" 4 Mpc" 1 and k = 10" 1 Mpc" 1 indicated 
on figure [5] correspond roughly to the quadrupole and the I = 2000 mode of the 
multipoles of the acoustic peaks. It is interesting to note that the first year WMAP 
data are equally consistent with both curves shown in figure |||; the solid one having a 
constant value of dn s /d In k, and the dashed one having dn s /dlnk = outside of the 
range —7.2 < ln(k ■ Mpc) < —3.8. Both of these spectra have the same likelihood, 
even without varying the cosmological parameters given in table [1]. We have verified 
this using our modified versions of CAMB and CosmoMC which allow for an arbitrary 
user-specified power spectrum. We will make further use this truncated version of 
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Figure 3: Log of best fit running spec- 
trum (solid, black) and the partial run- 
ning spectrum (dashed, red) versus In A;. 
The nonrunning spectrum is also shown 
(dot-dashed, green). 



Figure 4: Difference of x 2 between the 
nonrunning and running models, versus 
the minimum multipole included / m i n ; see 
eq. (2.3). The separate contributions 
from temperature (TT) and polarization 
(TE) as well as the total are shown. 



the running spectrum, which is flat at large and small k, referring to it as the partial 
running model. 

It is at first surprising that the partial running model should provide as good a 
fit to the data as the constant running model, since || showed that the evidence for 
running comes mainly from the lowest three multipoles. In contrast the running part 
of our partial running spectrum affects much higher multipoles. To better understand 
this, we have computed the difference in \ 2 between the nonrunning and running 
models including contributions from Z min to Z max (900 for TT and 450 for TE): 

Ax 2 (/ min ) = 5x 2 (nonrunning) — 5x 2 (running) (2.3) 

» — ^min 

where £x 2 (0 is the \ 2 contributed from diagonal (/) and off-diagonal (/' > I) terms. 
This is shown in figure f|. 3 We observe there that the low multipoles (I = 2,3,4) in 
the TT data do account for a substantial part of the total decrease in \ 2 between the 
nonrunning and running models; however the TE data, ignored by the analysis of |J, 
have the opposite effect and hence the A%f ot is rather insensitive to the inclusion of 
the first three multipoles. In fact, the larger part of the decrease in \ 2 clearly comes 
from the region near I = 40 in the TT data, which explains why our partial running 
spectrum starts running at relatively high k values, k = 0.001 Mpc -1 , compared 
to the the values which would affect the low multipoles, k = 10 -4 Mpc -1 . This 
perspective lends more interest to the possible confirmation or negation of large 

3 Since the original submission of this paper, we have reanalyzed the data using the WMAP 
three-year data release. The shape of the total A^ 2 curve is unchanged, but it has come down by 
approximately A\ 2 — 1 due to a similar shift in the contribution from TE. 
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running by future improvements in the data, since the experimental determination 
of the higher multipoles is not so limited by cosmic variance. 

Although the LSS data do not favor the running models as shown in table [| we 
have also looked for the best fit models using the partial running spectrum, where we 
allowed the low and high k values which define the partial running model to vary. 4 
Table |3| shows the best fit values and the likelihood for the different data sets. Since 
other CMB data strengthen the evidence for running, but LSS data weaken it, when 
including LSS data, the lower cutoff will be larger than, and the higher cutoff will 
be smaller than, those including other CMB data. All the x 2 values for the partial 
running model are seen to be lower than or equal to the running and nonrunning 
models. 
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Table 3: The likelihood of the partial running model, using the data described in table |2|. 



To summarize our view of the experimental evidence: there is an interesting hint 
of a running spectrum which is driven in large part by multipoles near I = 40 in 
the CMB data; the inclusion of large scale structure data do not reinforce this, but 
they also do not contradict it. Although the statistical significance is not great, it 
can also be argued that it is not negligible. According to the Bayesian information 
criterion [31], one should decide whether the addition of new parameters is justified 
by considering whether the combination 



BIC = x 2 + k\nN 



is decreased, where k is the number of parameters and N is the number of indepen- 
dent data points. Not all multipole moments of the CMB should be considered as 
independent, since it is known that the acoustic peaks can be fit with great accu- 
racy by a Fourier series with 20 terms. A change of \ 2 of 3 can therefore support 
3/ In 20 = 1 new parameters. Although we introduced two new parameters, the 
running and the tensor-to-scalar ratio, there is a degeneracy between these, which 
makes them more like one parameter. The effect we are discussing is therefore on 
the borderline of being significant, a situation that could change in either direction 
with new data. 

4 In trying to understand why the SDSS analysis marginally favors running, yet CosmoMC 
using SDSS data does not, we note that the SDSS analysis uses a non-flat model when running is 
considered. 
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3. Reconstruction of Potential 



3.1 Formalism 



To find inflaton potentials corresponding to a desired power spectrum, the Hamilton- 
Jacobi formulation of inflation is useful |23], ^4j. Here we introduce a method based on 
integration over H(k), which is simple and accurate; see ref. p2]| for more discussion 



and results in constant spectral index case. Ref. [21| presented a similar integration 
method, but used the lowest order slow-roll approximation in each step, and hence 
their method is not as accurate as Hamilton- Jacobi formulation which we describe. 
Initially, one takes the Hubble parameter to depend on the inflaton field, H = H(<p). 
In units where Mp\ = 1, the Friedmann equation can be written 



3H 2 - 2H 



12 



V 



where H' = dH/ dip. The inflaton velocity is given by 

6 = -2H' 



(3.1) 



(3.2) 



We begin by using the lowest order prediction for the power spectrum of the 
curvature (scalar field) perturbation in the slow-roll approximation: 



P(k) 



H 4 

c a 2 



(3.3) 



where C\ = 4-7T 2 , and the k dependence is implicitly defined using the horizon crossing 
condition, 

k/a = H — > Ink = lna + lnH (3.4) 

Since, using ( |3.2j ), In a = f* dtH = d(j)H/<j) = — | d(j) H/H' , one can take d/d(j) 
of eq. ( |3.4|) to obtain a first order differential equation determining k(<fi): 



dink H' H 



H 2H' 



(3.5) 



Our strategy for reconstructing the potential is to rewrite the above equations 
using In k as the independent variable instead of 0. By straightforward manipulation 
using the chain rule, and eliminating from eqs. ( |3.2| - [3.3|) , we find that 



dH 
dink 



dink 



H 3 



2c x P - H 2 
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1/2 



2Hy/^P 

~2 Cl P-H 2 



(3.6) 



and we choose the plus sign throughout this paper. This choice of sign has no 
physical significance during any period when <ft is changing monotonically with time, 
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as during inflation. These equations can be integrated numerically to find H and 
as functions of Ink. This implicitly defines H((p), which can then be used in (|3.1| ) to 
find V(4>). There is one relevant constant of integration H , the initial value of H, 
which determines the overall inflationary energy scale. The initial value of has no 
physical importance. For example we can take = when k is at the pivot point fco 
at which the power spectrum is normalized in ( |2.1| ). 

In the following, it is convenient to refer to rescaled variables. We define a 
rescaled power spectrum 

P = AV (3.7) 

where A is the amplitude of the spectrum in eq. ( |2.1| ), and thus V is the ex- 
ponential factor V = exp[(n s — l)ln^ + l^^hi 2 ^]- A is determined by the 
COBE normalization, which gives P(k c ) = 25/4 x (1.91 x 10~ 5 ) 2 ~ 2.3 x 10~ 9 
at k c = 7.5a H ~ 0.002MPC" 1 , consistent with the values in table [I]. We further 
rescale 

H = aH; a 2 = 2 Cl A (3.8) 
The equations of motion become 

n 3 



H 



V-H 2 



T v -re 

Similarly, we define the rescaled potential through 



(3.9) 



V = a 2 V = a 2 



o 2- 



n 



m 2 - 2{H/(j>) = a 2 H 2 ( 3 - ^- ) (3.10) 



The Hamilton- Jacobi slow-roll parameters are 



H ,s 2 



H 

From eqs. (^]2| - 13~^ ), we can rewrite the e as 

n 2 



2(-] (3.11) 

n = i[^r) (3.12) 



(3.13) 



Since e < 1 is required for accelerated expansion, Ti 2 is always smaller than V during 
inflation, so the equations of motion (|3.9|) remain nonsingular. 
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3.2 Reconstructed Potentials 

The method described above can be used to numerically determine a family of po- 
tentials V(4>) which give rise to a desired scalar power spectrum P(k). There is no 
unique potential for a given P(k) since changing the scale of inflation can be compen- 
sated by a change of the slope of the potential. This can be seen through the relation 
P ~ V 3 /V' 2 . Later on we will break this degeneracy by invoking information about 
the tensor power spectrum. 

As a first step, we reconstruct the potentials only over the region of field space 
traversed by the inflaton during the first ~ 8 e-foldings of inflation after horizon 
crossing, since this is the relevant period for affecting the CMB. We numerically 
integrated the equations ( |3.9| ) subject to the initial conditions H = Ho, 4> = 0. To 
reproduce the spectrum shown in figure p], we take the initial value of In fcj = —11.8, 
which corresponds to times several e-foldings before horizon crossing, and integrate 
to a final value oi Inkf = —1.4. The derived potentials, and the range over which 
<fi changes, depend strongly on the choice of Ho- However, fairly universal behavior 
can be seen by graphing the function 

m = ,3, 4) 

where 0f is the final value of 0, and V^f are respectively the initial and final values 
of V, which is a monotonically decreasing function. This is shown in figure |5|. The 
rather narrow range of potential shapes shown are actually inclusive, since potentials 
with smaller values of Ho are very close to the curve for Ho = 0.01, while the largest 
possible value of Ho (given by y/V(ki) = 0.215 using the spectrum (|2.1[)) yields a 
potential nearly coinciding with that of Ho = 0.2. Notice that Ho is not allowed to 
exceed Viki) in ( |3.9| ). This traces back to the physical restriction that the Hubble 
parameter can never increase during inflation, or equivalently, the kinetic energy of 
the scalar field cannot be negative. 

A much greater variation in the shapes of the rescaled potential can be observed 
if we insist on maintaining the running index power spectrum over a longer period 
of inflation, going beyond that which affects the observable multipoles of the CMB 
temperature fluctuations. Figure |6| shows that the relative shapes of inflaton potential 
in different initial conditions (lnfcj = —11.8, Inkf = 8.8). We can see that they are 
quite different, but all of them have a "bump." (Potentials with Ho smaller than 0.01 
are close to the curve of Ho = 0.01.) 

Although the relative shape of the inflaton potential is somewhat insensitive to 
Ho (at least over shorter periods of inflation), the scales over which it changes do 
depend strongly on Ho- The dependences 

<p { ~H ; V-^Hl Vi-Vf~^ (3.15) 
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Figure 5: Reconstructed potentials for 
a range of initial Hubble parameters. 
"Corrected" curves are discussed in sec- 
tion 3.3. 
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Figure 6: Reconstructed potentials of 
different initial Hubble parameters, inte- 
grating over 20 e-foldings of inflation. 
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Figure 7: Dependence of potential parameters fa, V;, and (V; — Vf), on the initial Hubble 
constant. Equations are numerical fits to the small-Wo parts of the curves. 



are shown in figure [?]. These dependences can be understood from equations (|3.1| - 
|373|). Since P(k) does not vary, eq. (|3.3|) implies the scaling H' ~ H 2 . Thus eq. ( |3.2| ) 
predicts that changes by an amount A<p ~ H 2 At ~ H, since we are interested in 
time scales of order 1/H. Eq. (|3~T| ) obviously predicts that the magnitude of V scales 
like H 2 . The variation in V is AV ~ ^'A0 ~ M'A0 ~ # 4 . 

The numerical value of 7io has no physical significance except in comparison to 
the value of \/V at the chosen starting time (i.e., starting value of Ink). Suppose we 
evolve the system for two different initial conditions, (lnk\, Tii) and (hi/^, 7"i2) with 
ki < k 2 . Given the solution for (lnfci, 7"ii), it is always possible to find (lnfc 2 , ^2) 
such that the two solutions match in the region where they overlap. However the 
converse is not true. This is because of the attractor nature of inflationary solutions: 
since the effects of variation of initial conditions are exponentially suppressed with 
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time, it is not generally possible to find initial conditions at very early times which 
match an arbitrary deviation from slow-roll behavior at later times. This point is 
important if inflation did not start much earlier than horizon crossing, so that the 
transient effects of initial conditions can be visible in the low-fc part of the power 
spectrum, as we will discuss in section |3T4\ 



3.3 Slow- Roll Corrections 

In the previous subsections we reconstructed the potential using the leading slow- 
roll approximation for the power spectrum. Here we show how corrections which are 
subleading in the slow- roll expansion can be incorporated in a simple way. 5 The next 
order result in the slow-roll approximation is p4| 



V {2 \k) = [l- e +(2-ln2- 7 )(2 e -77)] 2 P (0) =CV W (3.16) 

where 7 ~ 0.5772 is the Euler-Mascheroni constant, is the lowest order power 
spectrum which we already computed. To incorporate these corrections, which are 
dominated by 77, we make the replacement V — > V/C in eqs. ( |3.9| ), where the correc- 
tion factor C is computed using the previously obtained solution for Ti. This proce- 
dure can be iterated until the solution converges. Convergence to a self-consistent 
solution is obtained after three iterations. In this way we avoid having to solve a 

00 00 

higher order differential equation for Ti,. The second derivatives Ti, <f> in rj can be 



computed in terms of Ti, using the equations of motion (|3.9|) . The result is that 
the shapes (and magnitudes) of the potentials are slightly modified, as shown by the 
"corrected" curves in figure |5]. The maximum value of Tio is also reduced from 0.21 
to 0.11. The changes to the parameters in figure [F] are so slight that we do not show 
them. 

3.4 Sensitivity to Initial Conditions 

Once the potential has been determined, the power spectrum does not uniquely 
follow. One must also supply the initial velocity of the inflaton field. Normally this 
is unimportant since the inflaton reaches the attractor solution in which the slow- 
roll approximation 3H<j) = — V is satisfied, hence is determined by the potential. 
However the initial velocity required to match the desired power spectrum does not 
generally coincide with —V/3H; instead it is given by eq. ( |3.3j ). 

The effect is most pronounced for large values of Tio, as illustrated in figure |[ 
In this example, (where lnfcj = —9 and Tio = 0.55, and we do not include slow- 
roll corrections) the initial value of implied by eq. (|3.3j) is <ft = 0.76 (in Planck 



5 Recently ref. |33| investigated the accuracy of the slow-roll expansion for quantities like n s and 
dn s /dlnk in models with large running, finding that the inclusion of higher order corrections did 
not improve the fit of the approximate spectrum to the exact one. However, this leaves open the 



question as to whether the failure is with the approximation ( 3.16| ), or with its Taylor expansion in 
powers of In k. 
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units), while the slow-roll value is (f) = 0.33. Using the latter value, the spectrum 
is distorted in its shape at the beginning of inflation, and by a constant horizontal 
shift A In k thereafter. The figure shows that compensating for A In k causes the two 
spectra to match well once the transient effect of the initial condition has died away. 
For smaller inflationary scales (small TCo), the shift Aln/c becomes negligibly small, 
while the low-k distortion remains significant. Since the horizon-crossing condition 
is that Ink = In a + In H, this shift in In A; can be absorbed in a change of the scale 
of inflation. 




/ 

/ 

/ 

-1 / 
/ 
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Figure 8: Power spectra obtained from reconstructed potential, for different initial con- 
ditions of (f>. 



4. Tensor Contribution to CMB 

It was noted in section ^| that the power spectrum with large running provides a 
good fit to the data only if a large tensor component is also allowed. A degeneracy 
between these parameters is understandable, since large running suppresses power at 
low k, while the tensor contribution enhances it. In single-field inflation, which we 
are assuming, the tensor power is related to the scalar power, since at lowest order 
in the slow-roll approximation, 

P t (k) = —H 2 (4.1) 



(compare to the scalar power, eq. ( |3.3j )). This convention for the normalization of P t 
is consistent with the definition of the tensor-to-scalar ratio r used by CAMB and 
by ref. §, 

r = Wi = i6eih) - (4 - 2) 
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Parameterizing the scalar spectrum as eq. ( |2.1|) and using the WMAP parame- 
ters in table |], we can determine the tensor spectrum as well as the reconstructed 
potential uniquely. However, a problem arises with the running spectrum when we 
try to achieve a large tensor ratio r = 0.5 at the given pivot point k . Namely, eq. 
( |4.2|) shows that e ~ 1/32 at fco; however if we try to evolve backwards in time from 
this point using the running spectrum, the slow-roll parameters diverge, and prevent 
us from extending the initial value of In k to as small values as desired. This be- 
havior is illustrated in figure || which shows the reconstructed slow-roll parameters 
corresponding to the WMAP running spectrum. The figure shows that e remains 
small (< 0.05) from Ink ~ —7.2 to 1.3, rj runs linearly and gives the running of 
the power spectrum. But demanding a large tensor-to-scalar requires the slow-roll 
parameters to be quite large at the lower scales near k ~ 10 _4 Mpc _1 . Moreover, 
negative running at large k makes the power spectrum decrease exponentially, again 
leading to an increase in the slow-roll parameters, since e = 7i 2 /V. The slow- roll 
approximation breaks down before the end of inflation. For example, we get only 
15 e-foldings while e < 0.2 in the WMAP running model. These observations are 



consistent with ref. ||10|| , who noted that e reaches a minimum near the scale where 
n s crosses 1 in models with large running. 

To alleviate the above problem, and extend the duration of inflation before the 
pivot point, we use the partial running spectrum with a low-fc cutoff which was 
discussed in section 0, and shown (together with the constructed tensor spectrum) 
in figure |10[ We could also include the high-fc cutoff on the running, but we choose 
to omit it, in order to discuss the constraints which arise on the length of inflation 
if the spectrum continues to run at high k. The tensor spectrum is not a constant 
as is usually assumed when computing the CMB temperature anisotropy, and we 
have modified CosmoMC to use the constructed tensor spectrum. However, we do 
not observe any effect on the likelihood of the best fit running spectrum, nor any 
change in the other cosmological parameters, relative to the standard form of the 
tensor spectrum. 

With the low-fc cutoff on the running of the spectrum, the beginning of inflation 
can be pushed back before horizon crossing, while maintaining a large tensor ratio 
at the pivot point; this is shown in figure O, where the slow-roll parameters are seen 
to remain small during the relevant part of inflation, while the minimum value of e is 
consistent with r = 0.5. The discontinuity in 7] is due to the discontinuity in the slope 
of the power spectrum where it goes between nonrunning and running, which could 
be avoided by making a smoother transition. On the other hand, since we chose not 
to cut off the running at high k, there is a limit on the amount of inflation after 
horizon crossing. The reconstructed potential continues to decrease monotonically, 
so there is no abrupt end to inflation, but e becomes steadily larger, and quasi- 
exponential inflation crosses over to power-law inflation with a small power. This 
model gives 16 e-foldings of exponential inflation, during which e < e« = 0.12. The 
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Figure 9: The reconstructed slow-roll 
parameters for WMAP running power 
spectrum, showing singular behavior at 
low k. 
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Figure 10: The partial running scalar 
spectrum and the reconstructed tensor 
spectrum (unnormalized, A = 1). 



transition to reduced acceleration is shown in figure [12], where one can see a break in 
the logarithm of the scale factor versus time at around In a = 20 e-foldings. Figure 
13, which plots In a versus lnt, shows that after the break the expansion is indeed 
power-law. 
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Figure 11: The reconstructed slow-roll parameters for partial running power spectrum. 

With the given tensor ratio, we can uniquely determine the scale of inflation and 
the single out one reconstructed potential from the family that was derived in section 



372j . Members of this family were parametrized by the ratio 



«l = J^ = l^ = ^ = e (4 .3) 
V 2 Cl P 16 P 16 V ; 



evaluated at the initial value of Ink (not the pivot point). Figure |TT] shows that 



0.12 at the initial value lnfcj = —11.8, in order to get a large enough tensor 
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Figure 12: The logarithm of the scale 
factor as a function of time, showing the 
transition from exponential inflation to 
power-law expansion. 
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Figure 13: Same as figure 12, but plot- 
ting In a versus lni, to show the power- 
law character of the expansion at late 
times. 



ratio at the pivot point. We find that 7^o — 0.29, which satisfies the constraint 
7i 2 < V in the model with the partial running spectrum. We can then determine 
the scale of inflation; from table |], A = 2.1 x 10 -9 , we obtain a = \J1c\A ~ 
4.1 x 10- 4 M P1 ~ 9.9 x 10 14 GeV (see eq. (p|)). This corresponds to a Hubble scale 
of H = ali,Q = 2.9 x 10 14 GeV at the beginning of inflation, and an energy scale 
of V 1 ^ = (3H 2 Mffl A = 3.5 x 10 16 GeV. This is the energy scale expected in 
the slow-roll approximation with e ~ 0.1, V l/i ~ 0.027e 1 / 4 M P i ~ 3.7 x 10 16 GeV. 
The reconstructed Hubble parameter and potential, as functions of 0, are shown 
in figures YA and where we contrast the predictions of the leading order slow 
roll approximation and those including corrections at second order in the slow-roll 
parameters. The corrections are seen to be a small effect. 




Figure 14: The reconstructed Hubble 
parameter, at first and second order in 
the slow-roll approximation, for the par- 
tial running spectrum. 
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Figure 15: The reconstructed potential, 
at first and second order in the slow-roll 
approximation, for the partial running 
spectrum. 



The reconstructed potential shows the "bump" which is characteristic of models 
with large running from blue {n > 1) to red (n < 1), since a large third derivative of 
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the potential is required to get large running in the slow-roll prediction, 

dn„ 



(Jha : -2£v + 16e v w - 24e^ (4.4) 

in terms of the usual potential slow- roll parameters: ey = \(V' ' /V) 2 , r/y = V"/V, 
and £v — V'V" /V 2 . From the perspective of particle physics model-building, it is 
interesting to ask to what extent this potential can be matched by a renormalizable 
potential. We have therefore fit the reconstructed V to a quartic polynomial. Obvi- 
ously the fit becomes better as one restricts the range of the field, so we have done 
this exercise for a series of final values, which correspond to different maximum 
wave numbers in the spectrum. We made our fits using lnfcj = —9.2 (10~ 4 Mpc -1 ) to 
lnfcy = (1 Mpc -1 ), 3.8, 4.8, and 5.8, obtaining fits which for the case of hi fey = 4.8 
has the form 

V = 1.7 x 1(T 7 x (0.277 -0.163 + 0.0583 2 - 0.0102 3 + 0.00064 4 ) (4.5) 

where all coefficients (and 0) are in Planck units. Figure |16| plots the approximate 
fits against the actual reconstructed potential. Plotted in this way, all the fits look 
very good, but to better discriminate, we show the fractional deviations AV/V and 
correlation coefficients between the fits and the actual potential in figure [L7|. We 
see that final k values up to In kf = 4.8 provide good fits at the percent level. 
This In kf value corresponds to the first 18 e-foldings of inflation. The difficulty of 
sustaining 50 — 60 e-foldings of inflation with large running has been noted previously 
[§]. Our analysis gives further support to the idea that large running requires some 
significant change in the potential after the first few e-foldings, such as can occur in 
multifield models |12|, [13| . It is also noteworthy in our reconstruction that the inflaton 
changes by superPlanckian values during this initial period of inflation. However, the 
smallness of the coefficients in the potential indicate that the effective field theory 
description is not invalidated by the large field values. 



5. Conclusions 

We have revisited the question of a running spectral index, motivated by the WMAP 
first year data. We are allowing for the possibility that the large running applies only 
in a limited region of k space, where the CMB experiments are most sensitive and 
large scale structure is less sensitive. It is possible that even a parametrization with 
a running spectral index does not provide a good representation of the actual power 
spectrum over all observable scales. In fact we know that inflation cannot last for 60 
e-foldings if running remains at the level of dn/dlnk = —0.1 — we find only 15 — 20 
e-foldings of exponential inflation — so it is also theoretically well-motivated to fit the 
spectrum separately in low and high regions of k space. Something must change at 
larger k just so that inflation can continue for another 40 — 45 e-foldings. 
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Figure 16: The reconstructed inflaton 
potential, and a series of renormaliz- 
able approximate potentials, which are 
matched to the exact result over increas- 
ingly large ranges of field values (indi- 
cated by the wave number of horizon 
crossing corresponding to the maximum 
field value in each range). 
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Figure 17: The fractional errors of the 
fitting potentials. The correlation coeffi- 
cient r for each fit is also shown next to 
the label (r close to 1 corresponds to a 
good fit). 



Although the evidence for large running from the CMB is marginal, an im- 
provement in x 2 of 3, requiring the addition of two (although somewhat degenerate) 
parameters, we have argued that it is based on more than just the first few mul- 
tipoles of the temperature anisotropy; in fact, including the TE polarization data 
makes the first few multipoles unfavorable toward running. Rather, one observes a 
steady accumulation of A% 2 as more low multipoles are dropped, with a sharp drop 
occuring in the region of / near 40. Even though this is still in the region where 
errors in the WMAP data are limited by cosmic variance, it will be interesting to see 
if there any significant changes when the next two years of data will be released. 6 

Even though the large running scenario may be unlikely, it is interesting to pur- 
sue its consequences, given that it is not ruled out by the data. Since the slow-roll 
parameters do not remain very small over the entire inflationary trajectory, it may be 
useful to have a method of reconstruction of the inflaton potential which is not based 
on the slow- roll expansion around a particular pivot point, but instead satisfies the 
slow-roll prediction for the power spectrum (computed to any desired order in the 
slow-roll expansion) at every point in <fi or k space. We have devised an integration 
method of reconstructing V(<f>) directly from P(k), based on the Hamilton- Jacobi 
formalism, which has this property, and is also simple to implement numerically. We 
have shown how corrections to any order in the slow-roll expansion can be imple- 
mented without significantly complicating the algorithm. This method may be of 
interest in itself, independently of whether the running spectrum is confirmed by 

6 We are posting this paper on the astro-ph archive a day ahead of the second WMAP release, 
announced for 16 March 2006. 
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future data. 

We demonstrated the reconstruction method by deriving the family of inflaton 
potentials that reproduce the best-fit large running spectrum. We subsequently 
focused on the unique potential which gives a large enough tensor-to-scalar ratio r 
to be consistent with the best-fit spectrum. We pointed out that this is only possible 
when the running is cut off at low k values; otherwise r itself runs toward smaller 
values so quickly with k that it falls below the desired value r = 0.5 before the pivot 
point ko is reached. This is a new observation, though anticipated by [Hj . The large 
tensor contribution is one of the interesting features of the large running model; it 
puts the scale of inflation as high as possible (we find V 1 / 4 = 3.5 x 10 16 GeV) and it 
provides the greatest hope that direct evidence of tensor modes might be seen in the 
CMB. It also provides the best chance of uniquely fixing the inflaton potential, since 
without the tensor spectrum, and in the absence of a specific model, it is always 
possible to lower the scale of inflation while simultaneously making it flatter, to keep 
the scalar power fixed. We have shown that a renormalizable potential provides a 
good fit to the actual reconstructed potential over the limted region of field space 
where inflation occurs, even though the field values are Planckian. This provides 
another intriguing aspect to the large-running model: if true, it may give us an 
inflationary window on physics at the Planck scale. 
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